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Abstract 

The structure of critical lines of vanishing pion mass for the Wilson quark 
action is examined in quenched lattice QCD. Numerical evidence is presented that 
the critical lines spread into five branches beyond (5 = 5.6-5.7 at zero temperature. 
It is also shown that the critical lines disappear in the deconfined phase for the 
case of finite temperatures. 
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1 Introduction 



The problem of species doublers in lattice fermion formulations is well known. The 
Wilson quark action, commonly used in numerical simulations of lattice QCD, avoids 
this problem at the cost of explicit breaking of chiral symmetry. This means that the 
conventional Nambu-Goldstone mechanism associated with spontaneous breakdown of 
chiral symmetry cannot be invoked to predict the existence of massless pion for this 
quark action at a finite lattice spacing. It has been found, however, that massless pion 
does appear with the Wilson quark action once the bare quark mass mo, or equivalently 
the hopping parameter K = l/(2m + 8), is tuned appropriately. It is now widely 
accepted that the pion mass vanishes along a line K = K c (/3), called the critical line, 
which runs from K c (f3 = 0) ~ 1/4 in the strong-coupling limit (3 = Q/g 2 = to 
K c (/3 = oo) = 1/8 in the weak-coupling limit (3 = oo on the ((3, K) plane 0. 

One of us has suggested some years ago|2| that the massless pion which appears for 
the Wilson quark action may be understood as a zero mode of a second-order phase 
transition spontaneously breaking parity-flavor symmetry. It has also been argued that, 
while there are only two critical lines in the region of strong coupling, there exist ten 
critical lines in the weak-coupling region, paired into five branches which reach the weak- 
coupling limit /3 — oo at M = 1/2K = ±4, ±2 and 0. Some analytic[Q and numerical^] 
results supporting this view have been previously reported. 

More recently arguments have been presented that the end points of five branches 
move away from the weak-coupling limit for finite temporal lattice sizes corresponding 
to finite temperatures, so that the critical lines form five cusps at some (3 ^ oo, above 
which no critical line exists. The existence of the second critical line paired with the 
conventional one and a formation of a cusp by the two critical lines have been confirmed 
numerically in a recent full QCD simulation at finite temperature [||]. 

In this brief report we present further investigation of the structure of critical lines 
for the Wilson quark action using quenched QCD. One of the motivations of the present 
work is to find evidence for the existence of five pairs of critical lines toward weak 
coupling at zero temperature. We expect this structure to remain in the quenched 
approximation, in which case the structure should be computationally much easier to 
confirm than in full QCD. Another motivation is to examine how the structure of critical 
lines changes for a finite temporal lattice size N t corresponding to finite temperature. 
In particular we wish to examine how the first-order deconfinement phase transition, 
which takes place at a fixed value (3 = (3 c {N t ) independent of K for quenched QCD, 
affects the critical lines. Our expectation would be that massless pion, and hence also 
the critical lines, disappears in the high-temperature phase above the deconfinement 
transition (3 > f3 c {N t ). 
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2 Results at zero temperature 



Finding the location of critical lines throughout the (/3, K) plane requires a convenient 
indicator for their positions. For this purpose we employ the number of conjugate 
gradient iterations Nqg necessary to invert the Wilson fermion matrix to a given value 
of the residual. For infinite volume this quantity increases toward critical lines and is 
expected to diverge in between a pair of critical lines, where parity is spontaneously 
broken, due to the presence of zero modes of the Wilson fermion matrix. On a finite 
lattice where zero mode eigenvalues should be slightly shifted away from zero, we expect 
a peak in Nqq across each pair of critical lines. Let us note that the phase structure 
on the (/3, K) plane should be symmetric under K — ► —K. For K > we thus expect 
three peaks in Nqg as a function of K toward weak coupling. 

We have carried out a measurement of the number Nqq for a set of values of K > 
at (3 = 5.0, 5.5, 5.6, 5.7 on a 6 4 lattice, and at (3 — 5.8, 5.9, 6.0 on a 12 4 lattice. Results 
are obtained by an average over 60 (for 6 4 ) or 20 (for 12 4 ) configurations separated by 
1000 sweeps of the pseudo-heatbath algorithm. The stopping condition is taken to be 
r = y\\b — Dx\\ 2 /3V < 10 -6 where b represents a wall source of unit strength for all 
space-time sites and colors, and V = L 4 the lattice volume with L = 6 or 12. Runs 
are made only for K > since the structure of critical lines should be symmetric under 
K -> -K. 

Our results are shown in fig. 0(a) and (b). On a 6 4 lattice (fig. |l](a)) we observe 
appearance of three peaks at M m 3.0, 1.5 and starting at (3 — 5.6-5.7. On a 12 4 
lattice (fig. 0(b)) the three peaks become increasingly sharper as (3 increases. For the 
peaks at M ~ 3.0 and 1.5, the peak position shifts toward larger values of M, albeit 
only slightly over the range (3 = 5.5-6.0 examined in the simulation. These results are 
consistent with the expected structure of multiple critical lines appearing at (3 > 5.6-5.7. 

In order to confirm that the peaks of Nqq observed at (3 > 5.6-5.7 correspond 
to five critical lines of massless pion, we have calculated the pion mass at (3 = 6.0 
on a 10 3 x 20 and a 16 3 x 20 lattice. The pion propagator is averaged over 10-20 
configurations separated by 1000 pseudo-heatbath sweeps. The stopping condition for 
quark propagator calculation is taken to be r = \\D'b — D'Dx\\/\\x\\ < 10~ 10 with b 
a point source at the origin. We should note that we have encountered 'exceptional' 
configurations near the peaks of Nqq, on which the pion propagator takes a W shape as a 
function of time. Since such exceptional configurations are never encountered in our full 
QCD simulations^], we think that they are an artifact of the quenched approximation in 
which configurations with small fermionic eigenvalues are not suppressed contrary to full 
QCD. We have therefore excluded the exceptional configurations from our propagator 
average. 

We plot our results for the pion mass squared in fig. |2] as a function of M = 1/2K 
(see Table 1 for numerical values). A reasonable agreement of results for 10 3 and 16 3 
spatial sizes shows that finite spatial size effects are not severe in our data. We clearly 
observe that there exist four more critical values of vanishing pion mass in addition to the 
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Figure 1: Number of conjugate gradient iterations Nqg necessary to invert the Wilson 
fermion matrix with a wall source of unit strength for all sites and colors as a function 
of M — 1/2K. Stopping condition is chosen to be r = ^J\\b — Dx\\ 2 /3V < 10 -6 with 
V = L 4 the lattice volume, (a) 6 4 lattice for f3 = 5.0-5.7, (b) 12 4 lattice for f3 = 5.8-6.0. 
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Figure 2: Pion mass squared m\ at zero temperature as a function of M — 1/2K at 
(3 = 6.0. Open circles are for a 10 3 x 20 lattice, and filled circles for a 16 3 x 20 lattice. 
Lines show linear fits of ml in M. 
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conventional one located at the right most of the figure. Making a linear extrapolation 
of m% in 1/K using two or three data points, we find 1/2K C = 0.052(14), 1.577(63), 
1.653(5), 3.091(11) and 3.1848(44) for the five critical values from left to right. 

3 Results at finite temperature 

For a finite temporal lattice size N t corresponding to finite temperature, the critical 
line should be defined by vanishing of the pion screening mass extracted from the pion 
propagator for large spatial separations. In quenched QCD an interesting question is 
how the structure of critical lines changes across the deconfinement transition which 
takes place at /3 — (3 c {N t ) independent of the hopping parameter K of valence quark. In 
the confining phase (3 < /3 c (N t ), it is natural to expect that the structure of critical lines 
remains qualitatively the same as at zero temperature except for a possible shift of their 
location by a magnitude depending on N t . On the other hand, the pion screening mass 
would not vanish in the deconfined phase (3 > (3 c (N t ), and hence critical lines would 
disappear. For values of K close to the conventional critical value, such a behavior has 
been seen in a previous work0. Our aim here is to carry out a systematic study over 
the entire interval of the hopping parameter. 

Our simulation is made for two temporal lattice sizes N t = 4 and 8. The critical 
coupling is known to be (3 c (N t = 4) = 5.69226(41) §§ and f3 c (N t = 8) = 6.0609(9)0. 
For N t = 4, runs are made at (3 = 5.635 and 5.75 on a 12 3 x 4 lattice, and for N t = 8 at 
(3 = 5.90 and 6.10 on a 16 3 x 8 lattice. Quark propagators are computed on a periodically 
doubled lattice 12 2 x 24 x 4 in the former case and on the same lattice 16 3 x 8 for the 
latter. Numerical results for the pion screening mass are listed in Table 1. 

In fig. we plot by open circles the pion mass squared on a 12 2 x 24 x 4 lattice 
as a function of M — 1/2K at (3 = 5.635, which is in the confining phase. We find 
only a single critical point located at 1/2K C = 2.886(17). Exploring the region K > K c 
at K = 0.2-0.4, we found that the gauge configurations are dominated by exceptional 
ones for which the conjugate gradient solver for quark propagator takes over 5000-10000 
iterations to converge, while for K < K c a thousand iterations or less are sufficient. We 
take this as an indication that the parity-broken phase extends over oo > K > K c at 
j3 = 5.635 for N t = 4, i.e., the system is still in the strong coupling region where only 
one critical line exists for K > 0. 

For (3 = 5.75, which is in the deconfined phase, results for the pion screening mass 
squared are shown by filled circles in fig. |3|. While an overall pattern of m 2 as a function 
of M = 1/2K is similar to the zero-temperature case shown in fig. |2|, the pion mass does 
not vanish for any value of K for the present case. Thus the critical lines are absent in 
the deconfined phase as expected. 

Our results for Nt = 8 at (3 = 5.9 (confined phase) is plotted in fig. f|(a). We 
observe two critical values at 1/2K « 3.0. The pion mass also appears to vanish at 
1/2K « 1.65 — 1.66 and ~ 0.04, although our data is not sufficiently precise to resolve 
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Figure 3: Pion screening mass squared as a function of M — 1/2K on a 12 2 x 24 x 4 
lattice. Open circles are for (5 = 5.635, and filled circles for f3 = 5.750. Line shows a 
linear fit of in M. 
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Figure 4: Pion screening mass squared m\ as a function of M = 1/2K on a 16 3 x 8 
lattice, (a) (3 = 5.9, (b) f3 — 6.1. Lines in (a) show linear fits of ml in M. 



8 



if there are two critical values separated by a narrow gap or has a small but non- 
vanishing value. In fig. [|(b) we show how the behavior changes at /3 — 6.1 (deconfined 
phase). We clearly observe non- vanishing of pion mass in the deconfined phase. 

4 Conclusion 

In this brief report we have presented numerical evidence for the existence of five pairs 
of critical lines beyond (3 = 5.6-5.7 in quenched QCD at zero temperature, including an 
estimate of five values of K c at /3 — 6.0. 

It has recently been suggested that the two critical lines forming each pair merge at 
a finite value of /?, turning into a single linefllj. Our results show that this possibility 
is not realized up to (3 = 6.0 in quenched QCD. 

We have also shown that the critical lines are absent above the deconfinement phase 
transition (3 > j3 c (N t ). Since the deconfinement transition is of first order, we expect the 
critical lines to be sharply cutoff at the critical value (3 = j3 c (N t ). Our data, however, is 
not precise enough to examine this point in detail. 
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Table 1: Pion mass m T as a function of K. 
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